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LIST OF SYMBOLS 

vector used to relate p-^, pg* and v to the x,y,z 
coordinate reference, in. 
unit vector of A 
same as A, in. 
unit vector of B 

stall torque of servomotor for a given input voltage, oz-in. 

same as A, in. 

unit vector of C 

same as A, in. 

unit vector of D 

motor input voltage, V 

generator feedback voltage, V 

output voltage of pressure transducer signal conditioner, V 
error voltage of inner servo loop, V 

transfer function of amplifier-motor-generator servo loop 
transfer function of generalized control system with 
disturbances moved outside the loop 
ratio of servomotor torque to speed, oz-in- sec-rad 
unit vectors along coordinate axes of orthogonal reference 
systems (x,y,z) and (x^yS 2 ') 
imaginary quantity 
system inertia, oz-in-sec^-rad”-^ 
ratio of servomotor stall torque to input voltage, 



servoamplifier gain 

ratio of generator feedback voltage to shaft velocity., 
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V- sec-rad 

gear ratio of control angle (6) to motor shaft angle (0 m ) 

X r _ * 

ratio of the pressure difference between ports to -£he 
angular deviation of the port and velocity vectors from 
90 °, psi-rad“~ 

ratio of the output voltage of the pressure transducer 1 
signal conditioner to the pressure being sensed by the 
transducer, V-psi - 
system gain, sec”^ 

sensitivity factor relating angular deviation of port and 
velocity vectors from 90 ° to control angle changes 
system velocity-error constant, sec“^* 
torque delivered by the servomotor, oz-in 
sin 22.5° = 
cos 22.5° = 0 . 9211- 

pressure in pounds per square inch, differential 

pressure in pounds per square foot 
pressure difference between ports, psi 
reference input to servosystem, rad 
dynamic pressure, psi 
Laplaee-transform operator 
time variable, sec 
velocity unit vector 
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distance from (x,y,z) origin to the probe tip along the 
X-axis, in. 

pitch angle of probe (pitch up. = +, by definition) , rad 

* 

inverse of time constant due to pressure tubing lag, sec“~ 
angular motion used to position probe (control angle), rad 
positioning angle for null condition, rad 
angular deviation of 0 from its null position (i.e,, 

9 e = 0‘ - 0), rad 
motor shaft angular position, rad 

unit vector describing the orientation of the pressure ports 
system equivalent time constant, sec 

angular deviation of the port and velocity vectors from 
90°, rad 

portion of <f> due to the difference of the control angle 
from its null position (i.e., due- to 9 e ), rad 
portion of 0 due to a disturbance caused by motion of 
the other system's control angle, rad 
yaw angle of probe (yaw right = +, by definition), rad 
angular velocity of servomotor shaft, rad-sec""^ 
subscript 1 refers to all functions in the servo system 
controlling 9p 

subscript 2 refers to all functions in the servosystem 


controlling 02 



ABSTRACT 


Analyses and experimental results are presented on a control system 
used for positioning a mill seeking pressure prole . The probe is to 
be used in wind tunnels to measure wind velocity magnitude and 
direction. Two, interacting, positioning servomechanisms were used to 
control the pressure probe. The analysis of this system was accomplished 
using vector analysis as well as standard servoanalysis methods. These 
analyses showed the stability of the system to be dependent on the 
orientation of the null sensing pressure ports on the probe face. 
Laboratory and wind tunnel tests showed that the probe yaw and pitch 
angle measurements were accurate to within ±1° for angles up to ^5°. 

The development of this probe has resulted in a new useful instrument 
for wind tunnel research. 



CHAPTER I 


INTRODUCTION 

In wind tunnel research it has become more and more necessary to 
study the three-dimensional flow characteristics around the model. The 
magnitude and angular orientation of the airflow are of particular 
interest to the researchers. Various techniques and probe types have 
been designed to obtain this information (refs.' 1-|>) - In small tunnels, 
fixed probes are often used. Calibrated relations of the pressures 
between ports situated, at various places on the face of the probe yield 
the desired flow characteristics (refs . 6-10,). . In larger tunnels or' 
on aircraft automatic null-seeking probes can.be used- (ref. ( 11). Tbis 
type of -probe seeks out the direction of the wind velocity vector., and 
then by standard pitot techniques determines its magnitude. One of 
these probes has been developed at the National .Aeronautics and Space - , 
Administration’s Langley Research Center. 

It is necessary that measuring apparatus interrogating a point in 
space disturb the flow characteristics as little as possible . This 
specification along with those that follow presented a formidable 
problem- in the design of the probe. The other specifications are: 

1. The probe shall be capable of nulling on yaw and pitch angles 
up. to +4-5°. 

2. These angles shall be determined accurately to within ±1°'. 

3* Flow velocity magnitude shall be determined -accurately to within 

±1 per cent over a range of velocities corresponding to dynamic pressure 
between 2 and 200 psf . 


1 
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The probe tip shall maintain. a given point in space while 
searching for a null. 

This last specification made it necessary to design a unique positioning 
apparatus. A mechanism employing two angular motions which were not 
directly related to the yaw and pitch angles was designed to drive the 
probe to a given yaw and pitch configuration. The control of the, 
mechanism involves two, interacting, positioning servos, since one of 

> ' I 

the motions is superimposed on the other. .Most' conventional positioning 
servos are one -dimensional in nature, requiring just "one linear or ' 
angular motion. The control problem for this probe' is a little more 
complex because of the positioning interaction. Three-dimensional 
vector analysis techniques are needed to analyze the feedback control 
systems used. 

The purpose of this paper is to describe the techniques used in 
synthesizing the control systems. The mathematical analysis of the 
systems is also developed. Before discussing the control systems, a 
description is given of the probe and the method used for positioning 
it. This includes an analysis of how yaw and pitch angles are related 
to the angles of the controlling motions. The block diagrams and the . 
system equations are then developed. This is followed by the experimental 
results of wind tunnel and laboratory tests. A discussion of stability 
and error criteria is included. The paper concludes with a recommendation 
for future work and a commentary on the probe 1 s limitations and 


advantages. 



CHATTER II 


PROBE DESCRIPTION 

To ■understand the analysis of the control systems one must first - 
understand the method for positioning the' probe. There are five 
pressure ports on the hemispherical face of the pitot-type probe (refer 
to Figure l) . The axial port (5) is connected to one side of a 
differential pressure transducer, and the static ports on the body of. 

'the probe are connected to the other side of the transducer. This 

1 .. ( 

differential pressure is used, to yield flow velocity data after the 

probe has found a null and is aligned, with the' velocity vector.. - ’The 

other four ports consist of two sets of diagonal pairs ’( 1-2 and 3-t) . 

Each pair is connected to a differential* pressure, transducer which 

* * 

. > s 

provides error signals for one of the two positioning -servomechanisms 

. * • 1 ^ « 

* * 1 ^ , 9 ^ 

controlling ©^ and ©g. When the -prohe is in a nulled condition the 

. • * ■. 

pressures on the ports are equal; 4, and thus, the error' signals are zero. 
In this condition- axis 3 is colinear with the velocity vector. The 

< , Q 

axis of the pitot-type probe (axis 3) is the generatrix of a 22-1/2 
conical surface when revolved (0g) about axis 2. Axis 2 'is likewise, 
the generatrix of a second 22-l/2° conical surface when revolved ( 0p) 
about axis 1. This i-s the axis to which yaw and pitch angles are 
referenced. If the condition is met that the vertex angles of the two 
cones are equal, then rotations 9-^ and ©2 allow the probe to assume 
any orientation within the confines of a third cone. This cone has a 


3 
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45° vertex angle. Thus , the probe can assume yaw and pitch angles up 
±45° , while the probe tip maintains a given point in' space. 

Pressure lines run from the ports to the transducer’s located in 
the main body of the probe (.axis 1) . Flexible tubing is used to link 

' v 

' i » j ’ -* 

the fixed lines of the pressure ports to the transducer inputs {see 

1 * 1 

Figure 2).' Because of these pressure lines, the motion of 0^ must 

E *■ _ 4 • ■ 

be limited to jrl80°. It is also' necessary to limit the motion of 0p 
to ±l8o° due to electrical lines running' to the motor, generator and 
synchro used to control and sense the position of 0g. 

Yaw and pitch angle determinations . - .Figures 3 and. 4 illustrate 
the orthogonal coordinate systems which are used to relate, the control 
motions to the probe's yaw and pitch orientation (ref. 12). The 

AAA 

i - j - k coordinate system is fixed on the support axis of the probe 

AAA 

(axis l) . The 1’ - j ! - k' coordinate system is located a fixed 

AAA 

distance from the i - j - k origin as is shown by the vector A. 

The- orientation of the supporting member between axis 1 and axis 2 

A 

is represented by k'. Axis the axis of the sensing probe, is 

coline ar with C. The location of the tail of vector C with respect 

A A A - - 

to the i* - j ' - k* origin is identified by B. Axis 2 is, likewise-, 
represented by D and is colinear with i ' . The probe tip is situated 

* A A A 

a fixed distance x p from the origin of the i - j - k system. With 
the probe in the positions shown in Figure 3> are 

arbitrarily chosen to he 0°. Clockwise motions of 0p and &2 
have been chosen as positive, or increasing angles. Yaw and pitch 




I 




A in plane iVSc and 'i-k 
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angles are 0° for this orientation. Yaw right and pitch up of the 
sensing probe, with respect to axis 1, are considered as^ positive 
angle s . *■ , ' 

With = 0° and 0£ = 0° (see Figure 3) the , following relations 
can be obtained: 

■* 4 

- = - = sin 22,5° = 0.^84 = m (l). 

CD.- . ' 

*• « ¥ , 5 

- = - = ^ = cos 22-. 5° = 0.924 ='n (2) 

BCD- 

— * ‘ \ t, 

It is necessary to define the vector C, before yaw and pitch angle 
determinations can be made. When the probe has reached a null condition, 
this vector is colinear with the velocity vector. Referring to Figure 4, 
the relation between C, D, and- B is 

C = D - B. (5) 


and this, can he expressed as 


£ _ C _ D -B 
C C C 


or 



Since 


A 


a = 


i’ 


(*) 

(5) 


( 6 ) 
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and 


B = B sin 0o j ’ + B cos 0 2 k T 


(T) 


or 


h = sin 0 O j ' + cos 9o k' 


( 8 ) 


then 


c = n i' - m sin 0g j* - m cos 02 k'. 


(9) 


The unit vectors i', j', and k' must now he expressed in terms 

A A A 

of i., j ,, and k. Again referring to Figure it is seen that 


where 


1 *' = r— r = d 


D' 


( 10 ) 


D = Xp a -~ : k 


( 11 ) 


and 


A = - A -sin. 0^ j - A cos 0]_ k 


( 12 ) 


Then 


A Xp 1 + A sin 9 -l j + A' cos 0;^ k 
D 


A 

i ' - 


Using the relations in. equations (l) and (.2), 


(13) 


i 1 = n i + m sin 0-^ j + m cos 0-^ k 


(14) 
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If vector j 1 is perpendicular to the plane of i’ and A [i , 
k', and A are always coplanar), then 


A A 

9, i* x a 

0 = -=5S — 

i’ X a 


where 


a = — = - sin 8-, j - cos 0-, k' 
A ’ x -1- 


The cross product of i 1 and a’ can then "be expressed as 


, /N A 


1 0 


i’ X a = h m sin Bp ; m cos %l 


0 ‘ ‘ -sin 9p -cos 0p 


which yields 


i* x a = n cos 6p 0 - n/sin'0£, k 


The magnitude of the cross product is 


i* x a, = sin/ i 1 , a = sin(90° - 22.5°) = cos 22.5° = n (19) 


■Using the results of equations ( l8‘). and (19); 

/S A 

/\ n cos 0p 0 - n sin 0p k 


0 1 . = cos 0p 0 - -sin e-i_ k 


( 21 ) 
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The remaining direction vector k’ can "be expressed as 


i 1 ' x -j 


( 22 ). 


or 


k’ = 


,m sin 6^ m cos. 6-J 


0 cos 9-^ -sin. 0-]_ 


( 23 ) 


Solving the determinant, 


k J = - m i + n sin 8-, j + n cos 8-. k 


(2U) 


Using this equation along with (l4) , (2l) and (9), the unit vector 
c can he expressed as 

c = n^ i + mn sin 9^ j + mn cos 6^ k - m sin 9 0 cos G^q 

/s ^ A 

+ m sin 0g sin 0-^ k + cos 0 2 i - mn sin 9^ cos 0g j 

A. 

- mn cos 0.-j_ cos 0 2 k • (25) 

Collecting terms, 

l . \ ' * « 

c = (n^ + nr cos. 0 2 )i + (inn sin 0^ - m. cos G'^ sin 0 2 
- mn sin 9j_ cos 0 2 )j,+ ( 12 m cos 0^ + m sin 9^ sin 0 2 


- mn cos 9^ cos 0 2 )k 


(2 6 ) 


Since 


A «.* .A ' 1 A ’ ^ A> 

c = c x i : * c y j + c z k ' 


- (‘ 27 ) 
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then 


c x 


= n 2 + 


m cos Qc 


(28) 


Cy = mn sin 9^_ - m cos 0^ sin 9g - mri sin 0^ cos. 0g , ’ (2-9) 


and 

. ' ’ 'if 

c z - tan cos 0^_ + m sin 9]_ sin Sg' - mn cos 0]_ -cos 02 (30,) 

i 

t 

The yaw and pitch angles can he determined using equations (28-30), 
since c represents the orientation of the probe. The angle of attack 
(pitch) is represented by a, and can be computed from the’ expression 

c 

tan cl = — (pitch-up of the probe is positive) (31) 

c x 

The angle of sideslip (yaw) is represented by i}r, and can be computed 
from the expression 

-c 

y 

tan = - — (yaw right of the probe is positive). (32) 

x 

Figure 5 shows locus plots of probe yaw and pitch angles for 0-^ and 
02 values. 

It is to be noted that, since the velocity vector has a heading 
directly opposite- to that of the probe, the following relations exist: 

1. Yaw right of the probe is equivalent to yaw right of the 
velocity vector. 

2. Fitch up of the probe is equivalent to pitch down of the 


velocity vector. 



(+) Pitch angles 

('“) Pitch angles 

1 (+) Yaw angles 

— — - (-) Yaw angles 

Figure 5. Yaw and pitch angle plots. 








CHAPTER III 


DESCRIPTION OF CONTROL SYSTEM 

Manual and automatic control . --Figure 6 shows a schematic diagram 
of the "basic control system used for the two angular motions that 
position the sensing probe. The only differences between the two 
systems are gain values in the servoamplifiers, gear reduction ratios, 
and pressure sensitivity to probe motions. Pressures picked up at the 
probe face are transmitted to transducers in the body" of the probe . 

These transducers sense capacitance changes through motions of a stretched 
metal diaphram. The full range of the transducers is ±10 psid. They 
are so sensitive, however, that six other ranges are selectable by a 
switch on the signal conditioner. The most sensitive range is ±0.01 psid 
full scale. The signal conditioner provides amplification, so that a 
±5 Vdc output is available on any of the ranges selected. This output 
is chopped and fed to a servo amplifier , the push-pull output of which 
is the control voltage for a two-phase servomotor. An a-c tachometer 
driven in tandem by the motor- provides a damping signal back to the 
amplifier input. The motor drives its probe angular motion and a 
synchro transmitter through a high ratio gear reducer. The synchro 
provides a signal which is transmitted to a readout device used for 
monitoring and recording the angular position. This angular' readout 

% 

unit, located at the control chassis (see Figure 7), also- contains the 
limit switches whichprevent the probe from winding up and fouling its 
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Signal conditioners 
and power supply for 
pressure transducers 


Figure 7. Control panel. 
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pressure lines and control wires. Space and wire limitations prevented 
the limit switches .from being located within the probe itself. 

Component description . - Brief comments on the probe parts follow; 
A- ELECTRICAL COMPONENTS 


(l.) Pressure Transducers and Signal Conditioners - The pressure 
transducers used to' sense the pressures at the ports on the probe face 
are located in the body of the probe 1 s support axis ( see Figure 2) . 

These transducers employ stretched metal diaphrams., the motions of 
which are sensed as varying capacitance. This in turn modulates a 
10 kHz carrier voltage which is. supplied from a power supply located 
in the control panel. The signal conditioners (also located in the . 
control panel) provide the necessary gain and .demodulation circuits. 

Seven separate ranges or gain values may be selected by a switch on the 
signal conditioner. A ±5 Vdc output is provided for the full scale - ' 

k * l 

i ‘ 1 

value of each range. The selectable ranges are; +0.01 psid, ±0.05 psid, 

'i 

±0.1 psid", ±0-5 psid, ±1.0 psid, ±3-0 psid, an.d. ±10.0 ps'id. . i 


(2) ServoampTifiers - The amplifiers us.ed t to ' drive the control 
motors are solid state devices capable of delivering 5* 5 -watts from' a 
differential output. They are designed to drive the center tapped;/ 
control winding of a two-phase *K30 Hz servomotor. The gain can be set 
(up to 2500) by selection of an input resistor for single ended inputs. 

(5) Servomotors, Generators and Gearheads- - The control motor, 
feedback generator and gear reduction package are all integrated into 
one unit. Two units are provided and are identical with the exception 
of the gearheads. The 02 system has a 1*1-46:1 gear reduction, and the 



19 


9 -j_ system has a 1 - 96:1 gearhead with an additional !-:l external reduction. 
The motors operate with a 26 V, 1-00 Hz fixed phase and a center tapped 
control winding. Their no-load speed i-s 6500 xpm. The generators also 
require 26 V , 1-00 Hz and deliver 0.11 V rms/1000 rpm for equivalent 
damping. These damping signals are summed (l8o° out of phase) with the 
control phase signal at the servoamplif ier s , and the gains are set "by 
selection of the input resistors. 

(!■) Choppers (100 Hz) - Solid state choppers -driven at loo Hz 
convert the dc output voltages from the pressure transducer signal 
conditioners to signals compatible with the loo Hz servo components. 

(5) Synchro Transmitters - -The positions of the controlled 
angles are sensed hy loo Hz synchro -transmitters. The synchros are 
located in the prohe body, and are geared to the prohe drive axes-. 

(6) Angle Indicators - The, synchros in the prohe are wired hack 
to angle indicator units located in the- control panel. The synchro 
transmitter signals are fed into control transformers. These control 
transformers are driven to a null "by servomotors in the indicator units. 
Mechanical angle counters are geared to the servomotor -control trans- 
former systems. Limit switches are also geared to this mechanical 
system. These limit switches are wired to limit the travel of the 0-^ 
and 0p motions on the probe. 

(7) Control Chasis - The amplifiers, choppers, switches and 
relays used i-n controlling the servomotors are .mounted on a chassis- in 
the control panel (see Figure 7)* This unit also contains meters for 
monitoring the pressures on the null ports. 
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B« MECHANICAL COMPONENTS 

(1) Sensing Probe - ‘The pitot type sensing probe was constructed 
from 0.1875 in. o.d. stainless, steel tubing. Six smaller (-0.020 in. 
i.d.) stainless steel tubes are contained inside the large tube. The 
hemispherical face of the probe contains the ports for five of the 
smaller tubes. The total pressure port is located on the probe axis. 

The nulling ports are separated by 90° on a plane perpendicular to tbe 
probe axis. Each is situated on the probe face at an angle of ^5° to 
the probe axis. The static pressure ports are situated downstream 
from tbe probe face, and they feed, the- sixth small tube. 

(2) Probe Body - The .mounting assembly, the housings which 
contain the transducers, motors, generators, synchros and the drive 
trains were all designed at the National Aeronautics and .Space 
Administration’s Langley Research Center. The material used for the 
parts was stainless steel. 



TABLE I 


BIST OP COMPONENTS 


Component 

■Manufacturer 

Model number 

Pressure transducer 

Datame tries, Inc. 

- 1055 

Signal conditioner 

Datame tries, Inc. 

' 1015 

Power supply 

Datametrics, Inc. 

: 700 

400 -Hz chopper 

Solid State 
Electronics Corp. 

■ 6k • 

Servoamplif ie r 

Kearf ott 

C 7051 U 800 I ' 

Servomotor-generator- 

Clifton Precision 

73750-AAY-l 

gearhead 

Products 

arid 



. 73750-AAM-l 

Synchro transmitter 

Clifton Precision 
Products 

CGE-8-F-7’ . ‘ 

Angle indicator 

Servo Systems -Co. 

ID-1033 and. 
ID-103 1 ! 

Sensing Probe, Housings 
and Support Section 

HASA-LRC. 

— 
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Manual control .- The drive mode for the prohe motions is -selected 

by a switch located on the control chassis . When it is desired to 

* ‘ , > 

position the prohe manually, the input to the chopper is switched from 

t « * 

the transducer signal conditioner output to a battery and polarity, 

* V t ' 

switching network. The probe can be manually driven to a null by 
observing the metered outputs of the pressure transducers.' 

Automatic control.- As previously stated,' there are 'two angular 

■ — — — — — 1 ' i * \ . i ’ 

motions ( 0^ and 9 2 ) employed in positioning the probe to its proper 

5 t 

null. point. The motors controlling these motions are. driven by error 

r > * * 

signals obtained from differential pressure transducers and their i 
signal conditioners. One pressure transducer senses the pressure 
between ports 1 and 2 (Pp_ 2 ); und another senses the pressure between 
ports 3 and ^ (P^_4) • 

At first, it was thought that it would make- little difference 
which error signal was used to drive the motor controlling 9g. However, 
when 0 2 is .controlled hy the pressure difference signal at ports 1 
and 2,, the system is unstable. This can- he demonstrated by the 
following mathematical analysis: 

Referring to Figure 4, 

(1) The unit vector pp defines the orientation of a line 
which passes through point Xp and is parallel to a line between 
ports 1 and 2 (the direction of pp is from port 2 to port l) . 

(2) The unit vector p 2 defines the orientation of ports 3 

and 4 in a like manner (the direction of p 2 is from port 4 to port 3) • 
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The vectors and can he defined in terns of 9^ and 0^ 

as follows: 

From Figure ^ it can be seen that 


AAA 

p 2 = i’ X b * 


(550- 


Using equation (8)', 

P 2 = i 1 x : (sin 9 2 j ' +‘ cos k 1 ) 


m 

♦ 


or 


9 2 ~ - cos j ’* + sin 0 2 k' (55) 

Using equations (21)' and (25), 

A A A 

p 2 = - cos 0 2 (cos 9-^ j - sin 0-j_ k) 

+ sin 0 2 (- m i + n sin 0j_ j + n cos 0^ k) ( 56 )' 


and finally 

pg = - m- sin 0g i + (n sin 0^ sin 0 2 - cos 9^ cos 0 2 )j 

* /V 

+ . (n .cos 8]_ sin 0 2 + sin cos 9 2 )k (37) 

Vector p-j_ can be expressed as 

P x = p 2 x c ( 58 ) 


Using equations (55) and (9) 



2k 


A A 



Pi = 


0 


n 


-cos ©2 
-m sin Qq 


k* 

sin 02 
-m cos 02 


(39) 


Solving the determinant, 

A ^ q A A A 

P]_ = (m cos^ 0 2 + m sin^ 0 2 )i' + n sin e 2 ^ + n cos 0 2 k ’ C^) 


Using equations (l4), (21), and (24), 

A, A A A 

Pl = m(n i + m sin ©i j- + m cos 0i k) 

+ n sin Ogtcos j - sin k) 

4- n cos OgO- m- i + n sin ©i j + n cos ©i k) ’ (4l) 


Combining terms, 

,Pl = (inn - non cos 02-)i + (m2 sin 0^ + n cos 0i sin Qq 

+ n^ sin 0q cos 02)j + ( m ^ cos ®1 ** n S ^ EC ®1 s ?- n ®2 
+ n^ cos ©i cos 0g)k '(^2) . 

This defines p^, a unit vector which can he used to represent the 
orientation of ports 1 and 2. Likewise, P2 represents the orientation 
of ports 3 and k. The relationship of pi and pg to the velocity 
vector is now developed. 

.Let v = unit vector representing the orientation of wind velocity. 
For this analysis, assume that v is constant and only the prohe can 
move. At a null, 

v — — c (c = unit vector along the pitot prohe axis) (4j) 
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and 

v • p-j_ = v • p 2 = 0 (4^-) 

For small angular deviations, the, pressure differences between ports 1 
and 2 will be directly proportional to the amount that the angle 
between p^ and v differs from 90° • The same holds for the 
pressure difference 'between ports 5 and 4 and the angle between pp 
and v. This was shown by empirical testing of the hemispherical 
probe (see Appendix i),. The relation was found to be 

= 0.07 q.(psi)/deg (45) 

And q, <f>\} and are defined as follows: 

g_ = dynamic pressure in psi 

0-j_ = angular deviation of and v from orthogonality 

02 = angular deviation of p 2 and v from' orthogonality 

The dot products of v with p-^ and pg can.be expressed as 

v . p x = [vjlpjcos (9G° ± 0J) (46) 

and 

v • P 2 = |v| |p 2 | cos (90° t 0 2 ) (47) 

1*1 e N - N. " 1 ' w 



Since 
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these expressions reduce to 

v • p-j- = cos (90° t $]_) = + sin 0-j_ » + (^9)- 

and 


v . p g = cos (90° t 0 £ ) = * sin 0 2 « ¥ 0 2 (50) 

for small values of and 

■Determination of control motion sensitivities .- By determining 
B(v • p-j_) S(v * p 2 ) S^ 2 S(v • p-j_) 5(v • p 2 ) S0 2 

69 p d 9 p * • ^®2 d 9 p ^®2 ^®2 ^®1 ^®1 

(Refs. 13 and 14 ) the effect on pressure differences between ports 1 
and 2 , and ports 3 and 4 by Q^_ and 0 2 motions can be found. At 
any given null condition > v = - c, and 9-^ and 0 2 are at some fixed 
positions, 0 * and 9 ’ . The velocity vector can then be described 
from equation (26) as 


v = (- n^ - m^ cos 0 2 )i + (- mn sin '9£ + m cos 9-£ sin Q^’ 

A * ' * 

+ mn sin 9| cos 0^) j + (- mn cos 9^' - m sin 0£ sin 9 2 / * 

+ mn cos 9^ cos 0^)k , r \ *" • '51} 

* 

The effect of 0g on Pp _2 can found 'from 


i 4 t l 



+ 





(52) ’ 
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Since v is considered constant 


S _ 0 


and the equation simplifies to 


(53) 


acv • pp _ /^A + A Pj^ + 


30; 


’X 


to / na* 


(S0 2 




The terms v x# Vy, and v z can be obtained from equation (5l)> and 


B Plx 3p 


Pw, p 1v , and .p-, from (42). Then- after determining - 

' 302 302 


and 


30; 


the effect on Pq_g "by 9 2 some nul1 point can be 


expressed as 


3(v *• 0 L ) 


39/- 


9 1 - e i 

e 2 = 0’ 


= (- rr - tst cos 0£)(mn sin 9^) 

+ (- mn sin 9£ + m cos ; 0£ sin 0£ + mn sin 0£ cos 0£) 
(n cos 0£ cos 9£ - n 2 sin Q£ sin 8£) 

+ (- mri cos, '9^ m 'sin 9£ sin 8£ + .mri cos- 0£ cos -0£) 

1 i * 

(- n sin 9,£ cos 0£ - n 2 cos.0£ sin 9£) , , 

•' ? ( 55 ) 


Multiplying terms and simplifying ^ 

9(v • P X ) 


30 2 


= mn sin 0 * *cos 


[l -'(m 2 + n 2 )J t = 0 ( 56 ) 


Q 1 = e i 

0=0' 
2 2 



28 


' ■ APi o 

Since Pig is proportional to y • pi , then ~-is also zero. 

This means that Pp_g will not change for small changes in 9g from 

t)(v • pi )' 


its null position. By examining the function 


at values of 


©2=8^+ 8 £ 2 (®e2 = small angular change) 


e 2 ' 9 S - 0 S2 


3(» ■ P]j 

it can he seen that ^ = 0 at null conditions because a 

maximum exists-.' ;This results in a condition of instability. The -error 
signals due to must, be opposite in sign on each side of the null 

point.- This is not possible if a maximum exists. 


Since 9g cannot be controlled by pressures at ports 1 and 2 , 
B(v • pg) 

-~r — jg must be examined. This relation can be expressed 'as 


&(* • Pg) 


02x ft) + v *ft) + 




ft' 


rA 4. , r r p 2z 

P H^}, w 


And since 


,dv x _ = S ^z _ 

-Se 2 


( 58 ) 
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then 


d(v . p 2 ) /a P2x \ /5 P2y \ 

— 55T" = H^rj + v H 55 fj 


+ v z 



(59) 


Once again v x , v y; and v z are- obtained from equation (51) , but. 

^PPtt '^P?v 

P 2x’ ^2y’ an< l Pgg are found in (57) • After determining gg — , 

Sp2z 2 2 

and ■•, the effect on ^ by Qg motions at some null point 

can be expressed .as 


= (- n 2 - m 2 cos 9^)(- m cos 9p 
6- = 9| + (- mn sin 0^ + m cos 9^ sin 6^ + mn sin 0^ cos 0^) 

0 _ Q I 

2 2 (n sin 0^ cos 9^ + cos 9^ sin 9^) ■ 

+ (- mn cos 9^ - m sin sin 0^ + mn cos 9| cos 9^) 
(n cos 9£ cos 99, - sin 9J sin '9£) 


S(v . B 2 ) 
^2 


C6o) 


Multiplying terms and simplifying 


S(v . p 2 )’ 


'd9r 


= m 


*0-i = 


- 


e i 

9 2 


m 


This result shows that for smallnhanges in 9g about a’null position, 
Pj_ 4 _ will be directly proportional to changes in '02* It also shows 
that the sensitivity is independent of the null position ( 0 ^ 9^). 
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This analysis has shown that the motor driving 0 2 must be 
controlled by error signals due to pressure differences between ports 
3 and. 4. That leaves 0^_ to be controlled by pressure differences 
between ports 1 and 2. The sensitivity of Pp_2 1° ®1 motions can 

be determined from 


Since 


and 


B('v • P]_) - /Bv. 


Bg-i 


x\, / op lx\ 

= p lx(^7) + * 


\Be 3 r p iy\Be 3 J + v y(,Be 


y \ . r^xj' 

+ V 

"”“1 j 


l^zX ^ / Sp lz' 

+ Mad + T W, 


(62) 


Bv Bv Bv 
x _ y _ z 


BGn B9n B0n 


= 0 


C63) 


Bo 


lx 


Ben 


- = 0* since p-^ 'f f(0j-) 


m 


then 


B(v • pg_) ! /Bp 


~Bih 


= ,v- 


' v y®i 


lyi 


+ V, 


, fob*' 
HSeT. 


(65) 


The terms Vy and v z are .obtained' from equation (51) /-and pjy and 

, . Bpi^r ' ’ Bpn „ 

p lz from. (42). After finding and. , the effect on P^_ 2 

due to 0-j_ motions at some null point can he. expressed as 1 
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<Kv . 




e 2 = 


(- nm sin ,0£ + m cos 9£ sin 0^ + mn sin 0^ cos Q'J) 

(m 2 cos 0£ «• n sin ©■[ sin 0£, + n 2 'cos'0-[ cos >0£) 

+ (- mn cos 0£ - m sin 9J sin 9i, + Bin cos 9^ cos 0^) • 

(- m 2 sin 0-^ — n cos 0-[ sin 9£ - n 2 sin 0£. cos 0£) 

(66) 


Multiplying and simplifying, 

m sin 0^ (67) 

The following is revealed about changes in Pp_g for small changes in 
0-^ about a null position: 

(1) ‘ The magnitude and sign of the change are - dependent on the 
probe' 1 s null position due to the sin 0^ term. 

( 2 ) There are no pressure variations when 0 i, = 0 °, l8o°. 

It has been previously shown by equation (56) that the pressure 

between ports 1 and 2 is not disturbed by small changes in 0 2 about 

the null point. It will be of' interest to determine the disturbance 

on the pressure between ports 3 and ^ due to changes in 0 p. This 

c)(v . fe) 

determination is made by finding, — which can- be expressed as 

d(v . p 2 ) /Bv x \ 

= p2x W + 

+ p 2z(^j' + ‘ C 68 ) 



’+ p 



( bp 


2y\ Se 


+ v 


2y) 


y\ 60n 


^(v • p x ) 


60- 


0i = 

e„ = 


9 1 

9 2 
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The velocity vector (v) is constant, and thus. 


3v Sv Sv 

x _ y _ z 

§6^ 89^' S'0^ 


Since p-^ f f(0i)>. 


and equation (68) reduces to 


d(v * ,P2) VSp2y\ /^P2z^ 

— Se; y yl SeT + Msef", 


Terns Vy and v z are obtained from, equation (51) , and \ p 2 y and p 2 Z 

, . ^P2y ^P.2z 

from (37)* After finding r , and -V- — •> 1 

5e r . _de 1 


8(v • p 2 ) 
Sei 


= •(- mn sin +-m cos 0£. sin 0g,+ inn sin 0£ .cos 9g) 
'© 2 _ = &■[ (n cos 9£ sin 0£ ■fsin 0-^ cos 0g) 

Q =01 

2 2 + (- mn cos 0£ - m sin- 0£ sin 9g + mn cos 0£ cos 0i,) 

(- n sin 0^ sin 9£ + cos 9£ cos 0g) 


Multiplying and simplifying, 


8(v • p 2 ) 


= mn(l - 


■cos. e^) 
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The following is revealed about changes in for small changes in 

from, its null position-: ’ ► . 

(l) The change is dependent on the probe's null position due to 
the term (cos 91) 

'< 0 , • ) c 

(2.) The sensitivity i-s greatest at 91 = l80 . and 1 is aero- at 

' *■ > : * 

9 g = 0°.- ■ 

From this analysis it was found that there ’is a disturbance on the 
servosystem which controls Og due to the interaction caused b-y Qj, 
motions. However, the converse is not- true. Motions of 9^ (near the' 

s ,1 

null condition) do not affect the servosystem controlling 0j_. 



CHAPTER IV • 


ANALYSIS OF CONTROL SYSTEMS 


System equations and "block diagrams .- For any given 'velocity vector 
within the range of the prohe there is a set of yaw and pitch, angles 
which identifies its heading. These yaw and pitch angles are achieved 
through control motion settings of 0-^ and 0g. These angles (0£ and 
9g) are the reference inputs te the two feedback control systems. A 
block diagram showing the two control systems and their interactions 
for small signal analysis is shown, in Figure 8. 

The blocks G|(S) and GgCs) represent the cascade transfer 
.functions of inner loops. Each of these loops is made up of an amplifier, 
a motor and a feedback generator. The block diagram for this inner loop' 

is shown in Figure 9- The development of the transform G ' ( S) is given 

* » 

in Appendix II (Ref. 15). _ • 

4 * 

Each of the systems can he considered a simple positioning .servo- 
mechanism hy combining the cascade blocks and moving the disturbance 

I , ■ r 

outside the loop (ref. l6). Figure 10. is a generalized block' diagram 

* i 

of such a system. A new reference (R')' i- s formed by the effect of the 

. ,* ■ * * 

disturbance input. The transfer function for this system would be: ; ’ 




9(S.). 

K'(S) 


% 


S3 + 


h + kk„k, 


afg 


(S + 7) + Kj, 


(74) 
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Figure 8. Block diagram of control system. 
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where , 


Krji — 


WgaV 

J 


( 75 ) 


Discussion of the control systems .- The forward loop gain as taken 
from Figure 10 is 


G(S) =■ 


% 


sfs + il^aVs + 7) 


( 76 ) 


This can be expressed in terms of the velocity constant Ky and the 
equivalent time constant r as 


where 


G(S) 


V 


S(S + ij(S + 7) 


( 77 ) 


% - 


Kj 


h + kk a k g ) 


( 78 ) 


and 


T 


' J 

h + kk £ kg 


( 79 ) 


Since the system has unity feedback, G(s)' is. also the open-loop 
transfer function. The single pole at the origin makes this a type 1 

„ l • 

system. Theoretically, type 1 systems have -no steady state errors for 
a step input. This is desirable in a null seeking system such as this.. 

t ' ) ' * 
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It was not necessary to design for a low velocity error, because the 
probe is to be used for making static interrogations of points in a wind 
tunnel. In using the probe, turbulent conditions will often be 
encountered. To filter out the effects of these noise signals, the 
system’s variable constants can be adjusted to make the system a low 
pas's filter. Also a large damping factor is. necessary, since the angle 
indicator unit has a limited slewing speed. If turbulence is not a 
problem, and if a faster responding angle indicator were used (all 
electronic), the response of the system could be improved. Increasing 
the gain, using phase lead compensation, or a combination of both, would 
accomplish this. There is considerable gain margin- in each of the 
systems. This is illustrated in the root locus plots which are developed 
in Appendix III. 

9p system.- From Figure 8 it is obvious that while the 0p ■ system 

has no direct disturbances from- the Go 'system at a- null condition, 

t . , V 

it is affected by the ,0g system. This is due .to the fact that the 

. 4 

term sin 9’ appears in the k._ . transfer .function.' Two 'problems 
^ V 1 , * ^ * 

r * * * 

are immediately evident. First, a condition of' instability will occur 
if 9g is allowed to take on negative, values-. Secondly, the system 

f o , * * 

becomes an open loop when Qg =‘0- . The first problem was. solved by 

' \ 1 ' < 

limiting 9^ to positive values by setting, the limit switches at-: 

0 2 = 0° and 0g = 178°. This does not compromise the probe's range of 

' \ 

yaw and pitch angles. As can he seen in Figure 5, there are two null 
points for each combination of yaw and pitch angles within the probe’s 
range. One null point is in the negative 9g range, and one is in the 
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positive range . The second problem occurs only when 9 2 = 0 . It is 
true that 9^ is subject to random drifting 'when this occurs. However, 

when ■ 9 2 = 0°, the special condition exists that yaw ’and pitch angles 

* 

are both equal to zero. In this special condition 9^ can assume any 
position without loss of accuracy. This can be seen by again examining 
Figure 5. 

' - - ' * 1 

9 2 system.- This, system has no position dependent transfer functions 
- ' T " ^ • ' , > * ^ 
within the closed loop. However, it is subject to .disturbances from 

. k k * , ’ 

the _ 9 1 system. These disturbance signals, are dependent on the* position 

f ; ‘ . /' 

of 9 2 due to the term cos 9^ , in the .disturbance transfer function 
( see Figure 8j'. It is most fortunate that when - 9^ •, tends to drift 

i * V * 

(92 ='0°), the transfer function for its disturbing effects on 02 
takes on a value of zero- 



CHAPTER V 


EXPERIMENTAL RESULTS AND CONCLUSIONS 

Wind tunnel and" lab tests .- Calibration tests were made on the 
probe at Langley's 7 x 10-foot research facility model tunnel.. These' 
tests showed the probe- to have yaw and pitch angle errors of up to 2. 5°., 
The repeatability of the probe, however, was better than 0.2°. The 
errors seemed to be systematic, so further tests Were made with the 
probe in the laboratory. The results of these tests showed that the 
errors were due to two things: 

(l)- Deadband in the 0p drive gears of nearly 2°. 

(•2) Misalignment of the nulling pbrts. 

Both pp and p2 have apparent deviations of 6. 75° from the ideal 
condition of orthogonality with the probe axis. The word "apparent" 
is used, since the condition could be ’ caused by other factors such as 
an uneven pressure profile on the probe's face. 

Recommendations and concluding remarks ..- The deadband that was . 
found in the Gp system was removed by making some mechanical improve- 
ments on a -gear key.. The misalignment of the 1 ‘nulling ports .can be 
compensated for by setting an offset voltage on the signal conditioners 
for the null pressure transducers. The level of the offset will be 
dependent on the tunnel dynamic pressure (q.). The probe is presently 
waiting access to the model tunnel for further calibration tests. Even 
if the above corrections do not bring the absolute error to within 1° 
it is felt that the design objectives will have been met. The excellent 


4o 
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repeatability as shown in the first calibration tests insure that the 
absolute errors can be calibrated out. The sensitivities of the systems 
are such that the design requirement of a q range of 2 to 5 200- .psf was 
easily met and surpassed. The pressure ^-transducer sensing dynamic 
pressure by standard pitot tube techniques has an' accuracy of -better 
than- 1 per cent. This insures a determination of • velocity magnitude % 

1 * 

that is within the pi~escribed jri per cent. It is difficult to 'test the 

? , I > J t 

probe fully without the use of a tunnel. Plans are to use the probe 

in Langley's new V-STOL wind tunnel when that facility becomes 

* . * 

operational. The probe ''s full capabilities can best be examined at 
that time. 

If a similar probe were to be constructed in the future , several 
recommendations should be -offered. First, consideration should be given 
to the use of slip rings for the wires going to servo components. 

With slip rings, motions of Op need not be limited. Finally, an 
investigation could be made into the possibility of a more ideal 
orientation of the nulling ports on the probe face. Sensitivities 
could be improved, or cross coupled disturbances might be reduced. 

In conclusion, it is felt that a useful instrument has been 
developed. While it is too large to be used in small wind tunnels, it 
should find application in many of the larger tunnels where subsonic 


investigations are made. 
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APPENDIX I 


PRESSURE SENSITIVITY OF NULL SENSING PORTS 

Tests were made in a model wind-tunnel to determine the 
sensitivity of the null ports on the prohe face. The probe was held 
fixed in a sting which could be accurately positioned in angle of attack. 
Port holes 1 and 2 were checked, and then the probe was rotated 90° for 
tests on holes 3 and L. The sting was driven through offset angles 
(j f ) ) of ±U°. Figures II -1 and II-2 are plots of differential pressures 
between ports versus offset angles. 

From these plots the sensitivity of both pairs of holes can he 

seen to he, 1 x 10 -2 psi per degree offset. Since the tests were made at a 

-2 

q. of 13.9 X 10 psi, the sensitivity can he expressed as: 




X 10 


psi/deg)( 


,13-9 X 10" 2 



or 


P p 

■ °- 07 

Tests were run at other q. values, and the hole sensitivity was always 
found to be 0.07 q(psi)/deg or ^ q(psi)/rad. 






APPENDIX II 


DEVELOPMENT OF THE INNER 'LOOP TRANSFER FUNCTION • 

The speed-torque characteristics of a two-phase ■ servomotor can ■ 
.be approximated by linear curves. A graphical representation of their 
relationship is shown in Figure II-l. The mathematical relationship is 

L = - h% + b 0 (II -1) 

.Over a limited range the ratio of stall torque to applied voltage can 
also he considered linear. Figure II-2 illustrates this. , From this 
curve 


Lq = ke a 


(n-2) 


is obtained. From Figure II-l and' equation (.II-l) it can be seen that' 


^o ~ ^o 


(H-3) 


and therefore. 


L = - + ke £ 


If the load is purely inertial 


L = J 


d 2 e 

dt^ 


:m 


and 


d 2 0. 


m 


d0 m 

= - h — St + ke 


at‘ 


dt 


k6 


(Il-h) 


(n-5) 


(II-6) 



hq, 



Speed (U) m ) 

Figure II- 1. Servomotor speed- torque curves. 



Figure II-2. Stall torque versus input voltage. 
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The LaPlace form of this equation is 

js 2 e m (s) = ke a (s) - hse m (s) 

or 

(JS 2 + hS)e m (S) - ke a (s). 

The transform of the motor can then he given as 

e mCs) k 

e a( s ) JS 2 + bS 

' - I { 

Using the relations' ‘obtained from Figure 9 in the text,* 

e a( s ) = k a e 6 (s) 
e e (S) = e t (.S) - e g (s) 

and, 

e g (S) = k g S0 m (S) 

The inner loop transfer function is 


(H-7) 


i (II-8) 


(iI-9) 


(xi-io) 

r 

, ( 11 - 11 ). 


( 11 - 12 )' 


G'(S) 


= 9 m( s ) 
e t( s ) 


J 


h + j 


S S + 


(II-13) 



APPENDIX III 


ROOT LOCUS PLOTS AND SYSTEM CONSTANTS 

The open loop transfer function , Ky and t -were described, 
in the text (equations 77 - 79) as > 


and 


V 


gCs)h(s) = g(s) =■ 


S S + i)(s + 7) 


(III-l) 


% = 


Kip , 


h + ^a k gl 


(III-2) 


where 


% 7 


(ill -5) 


and 


T = 


h + Kk k 
a g 


(III-4 ) 


The closed loop transfer function given in the canonical form., is 


G 


T 


R 1 + GH 2 

S 5 + ( 7 + —IS 1 " + J- S + 

1 r T • t ! 


iV + 1 a + 5T 


The characteristic equation is 


p + ^ + i^g 2 + 1 s + ~~ - 0 


' (III-5) 


(III- 6 )' 


L9 
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The closed loop poles are then found from the following constants: 

8 1 System 

The 0p system constants are: 

= 0.J84 sin 9^ 
k p L = 4 q. psi/rad 

kp = 5 /q . volt/psi (minimum transducer gain for 

prohe threshold at 0p = 0.5°) 

kp = 10"" 2 oz-in/volt 
* 

N. = 6k 



= 5-o4 x 10“^ 


Jl - 

: 8.35 x 10"^ oz-in- 

sec 2 /rad 


= 1.47 X. 10 ^ sin 0 ^ 

sec “5 

h l = 

: 4.0 x 10" il ' oz-in~sec /rad 

* 

k si 

- 3 x 10" 5 volt- sec 

/rad 

K ?1 

= I.06 sin 8^ sec”"' - 


T 1 = 

3-S x 10" 5 sec 


?1 = 

50 sec -1 



*These constants are adjustable. 
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The open loop poles are 

■S = 0, S = -50, and S = -278 

/ 

The closed loop poles are 

s = 0, -50, and -278 (0£ = 0) 

S = -1.1, -48-7, and -278-7 (e£ = 90°) 

The root locus plot for the 0 i system with 0£ = 90° is shown in 
Figure III-l. As 0g approaches 0° or' l8o° from 90° the closed loop 
poles approach the open loop poles. 

9g System 

The 0p- system constants are: 

&v 2 = 0-58^ 

k p 2 = ^ q. psi/rad 

k* =5/o volt/psi (minimum transducer gain for 

prohe threshold at 0 2 = 0-5°) 

k 2 = 10~ 2 oz-i-n/volt 

kf = 25 
a 2 

k =6-91x10-^ 
n 2 

-6 

J 2 =8.55 X 10 - oz-in-sec/rad 
Kj 2 = 7-9 x 10^ sec“5 
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hg - 4.0 X 10“^ oz-in-sec/rad 
kg 2 = 4 x 10-5 volt- sec/rad. 

Ky 2 =0.94 sec -1 
t 2 = 5*96 x 10"5 sec 

7g - 50 sec - "^ 

*These constants are adjustable. 

The open loop poles are 

S = 0, -50, and- -168 

The closed loop poles are 

S = - 0 . 96 ^ -48, 6, and -168.4 

The root locus plot for the 9g system is shown in Figure III-2. 







- Normalized 


FIGURE 3 
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APOLLO 9 SC104/LM3-RPS FILTERED DATA 
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foldout frame 












redout frame! I 


-gJLDOUT FRA^ 
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366610.07 366830.07 366850.07 366870.07 366860-07 - 366310.07 366930.07 36S9S0.07 366970.07 

TIME C SECONDS) 


366960.07 i 367010.07 



Mixture Ratio 



.49 











PRESSURE,,' REGULATOR' OUT (RSIA) 



TIME, SEC. 



Chamber pressure, psia 





Chamber pressure, psia 



Time, hr:min;sec 





APPENDIX A 


FLIGHT DATA 

Figure 

A-l APS Thrust Chamber Pressure (GP2Q10P - pCM) During BTD 

A-2 APS Oxidizer Isolation Valve Inlet Pressure (GP1503P-PCM) During BTD 

A-3 Fuel Isolation Valve Inlet Pressure (GP1501P-FM) During BTD 

A-4 APS Fuel Isolation Valve Inlet Pressure (GP1501P-PCM) During BTD 

A- 5 APS Fuel Tank Bulk Temperature (GP0718T-PCM) During BTD 

A-6 APS Oxidizer Tank Bulk Temperature (GP1218T-PCM) During BTD 

A-7 APS Helium Supply Tank No. 2 Temperature (GP0202T-PCM) During BTD 

A- 8 APS Helium. Supply Tank No. 1 Temperature (GP0201T-PCM) During BTD 

A-9 APS Regulator Out Manifold Pressure (GP0025P-PCM) During BTD 

A-10 APS Regulator Out Manifold Pressure (GP0018P-PCM) During BTD 

A- 11 APS Helium Supply Tank No. 2 Pressure (GP0002P-PCM) During BTD 

A-l 2 APS Helium Supply Tank No. 1 Pressure (GP0001 P— PCM) During BTD 

A-l 3 APS Fuel Differential Pressure - Tank Bottom to Engine Interface 
(GP061 6P-PAM) During BTD 

A- 14 APS Oxidizer Differential Pressure - Tank Bottom to Engine Interface 
(GP1116P-PAM) During BTD 

A- 1 5 Regulator Out Manifold Pressure (GP0025P-FM) During BTD 
A— 1 6 Oxidizer Isolation Valve Inlet Pressure (GP1503P-FM) During BTD 
A— 1 7 Oxidizer Injector Inlet Pressure (GP2006P-FM) During BTD 
A-18 Oxidizer Isolation Valve Inlet Pressure (GP1503P-FM) During BTD 
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